arXiv:1507.03756vl [nucl-th] 14Jul2015 


Relativistic Effects in 3-Nucleon Forces for Nuclear Matter and Finite Nuclei 


A.H. Lippok and H. Miither 
Institut fiir Theoretische Physik, 

Universitdt Tubingen, D-72076 Tubingen, Germany 

In order to simulate the relativistic effects of the Dirac Brueckner Hartree Fock approach for finite 
nuclei the part of the Urbana 3 nucleon (3N) force is considered, which represents the enhancement 
of the small components of the Dirac spinors for the nucleons in the nuclear medium. This 3N 
force is included in a Brueckner Hartree Fock calculation with rearrangement terms using a realistic 
model for the NN interaction. The strength of the 3N force is adjusted to reproduce the empirical 
saturation point of nuclear matter and then used in corresponding studies of the closed shell nuclei 
and "^^Ca. Special attention is paid to a consistent treatment of the spectrum of particle states 
in the NN propagator of the Bethe-Goldstone equation. 

PACS numbers: 21.60.Jz, 21.30.Fe, 21.65.-f, 26.60.Gj 


I. INTRODUCTION 


One of the main challenges of theoretical nuclear 
physics is the attempt to derive the bulk properties of 
nuclear systems, which includes the saturation proper¬ 
ties of infinite nuclear matter but also the binding en¬ 
ergies and size of finite nuclei, from a realistic model of 
the nucleon-nucleon (NN) interaction. In this context 
models of the NN interactions are called realistic if they 
have been determined to fit the experimental data for 
2 nucleons, i.e. the NN scattering phase shifts and the 
data of the deuteron, with high accuracy. These so-called 
realistic NN interactions are in contrast to phenomeno- 
logkal or effective NN forces like the S^rme interactions 
[ll,0 or relativistic mean field models Q, which are fitted 
to describe the bulk properties of nuclear systems using 
Hartree-Fock or mean-field approximations. 

Examples for such realistic NN interactions are the 
One-Boson-Exchange (OBE) Models of the Bonn (re¬ 
spectively Idaho) 0 group or the local interaction mod¬ 
els of the Argonne or Urbana groups Rather so¬ 
phisticated versions of these models have been devel¬ 
oped like the CDBONN potential^ or the Argonne V18 
potential^. 

Although these potentials are “soft” as compared to 
the hard-core potentials developed in the middle of the 
last century, they contain strong tensor- and short-range 
components, which make it inevitable to employ a calcu¬ 
lation tool which treats correlations beyond the mean- 
field or Hartree-Fock approximation^. A typical ex¬ 
ample of such a many-body approach is the Brueckner- 
Hartree-Fock (BHF) approximation, which is based on 
the solution of two-nucleon scattering equation in the 
nuclear medium, leading to an energy- and density de¬ 
pendent effective interaction, the so-called G-matrix. 

Attempts have been made to compare this G-matrix 
with the phenomenological models of a nuclear force as 
the Skyrme interaction mentioned above and identify the 
density dependence of the Skyrme interaction with the 
medium dependence of the G-matrix[l, Q. 

The variation principle of Hartree-Fock calculations 


with density-dependent forces leads to rearrangement 
terms, which in the case of the Skyrme interaction are 
very important to obtain good agreement with the em¬ 
pirical data for nuclear matter and finite nuclei. There¬ 
fore also the BHF approximation has been extended to 
include such rearrangement terms. Accounting for the 
energy dependence of G this leads to the Renormal¬ 
ized BHF (RBHF) approach [l^ whereas the Density de¬ 
pendent HF (DHF) approximation also accounts for the 
Pauli-rearrangement terms [^, [ll|. The inclusion of the 
rearrangement terms is not only justified to simulate the 
features contained in the effective theory, but, in contrast 
to BHF, the DHF approach fulfills the Hughenholtz van 
Hove theorem due to the inclusion of the rearrangement 
terms. 


m 


The BHF approach using realistic NN interaction leads 
to a saturation point for nuclear matter. The saturation 
points calculated with various models for the realistic NN 
interaction form the so-called Goester band(^.fl3l|. which 
misses the empirical data in a significant way. There¬ 
fore many attempts have been made to explore various 
effects, which make the NN interaction in the nuclear 
medium different from the NN interaction in the vac¬ 
uum. An example of such studies is the inclusion of sub- 
nucleonic degrees of freedom e.g. in terms of intermedi¬ 
ate excitations of the interacting nucleons to the A(3,3) 
resonance. Such mutual polarization effects should occur 
already in the interaction of 2 nucleons in the vacuum. In 
fact such processes with intermediate isobar excitations 
provide a substantial contribution to the medium range 
attraction of the NN interaction, which in the OBE model 
is described in terms of the exchange of a light scalar a 
meson [Hill- In the nuclear medium these attractive 
terms are quenched due to dispersive and Pauli effects, 
features which could be described in terms of a density- 
dependent NN force or a 3-N interaction[I^ [I^. Such 
isobar effects turned out to be non-negligible, however, 
did not cure the problem of the Coester band, i.e. did not 
shift the saturation point calculated for nuclear matter 
to the empirical data. 

Another form of medium-dependence of the NN in¬ 
teraction has been supplied by the Dirac-Brueckner- 
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Hartree-Fock (DBHF) approximation (Tsl [T9l|. Here one 
considers the Dirac structure of the self-energy of the nu¬ 
cleon and accounts for the effect that the attractive scalar 
component of this self-energy yields Dirac spinors for the 
nucleons in the nuclear medium with an enhanced small 
component as compared to the corresponding spinor for 
a nucleon in the vacuum. Again this change of the Dirac 
spinors in the medium and the corresponding change for 
the matrix elements of meson exchange can either be de¬ 
scribed in terms of a medium-dependent NN force or by 
means of a 3N force parameterized in form of a Z-diagram 
with 2 (T exchange terms [2^ [2l| . 

The DBHF approximation has been successful in the 
sense that it reproduced the empirical saturation point 
of nuclear matter without adjustment of any additional 
parameter (l^[2l - [^ . Since however the consistent treat¬ 
ment of correlation and relativistic effects for finite sys¬ 
tems is a rather involved problem, the full Dirac Brueck- 
ner equations have not yet been solved for finite nuclei. 
In fact, Van Giai et al. addressed this problem as one 
of the main open problems in Nuclear Physics. Different 
approximation schemes have been developed, which treat 
either the relativistic effects or the correlation effects in 
an approximate way. 

The most popular approximation scheme is to ana¬ 
lyze the DBHF calculations of infinite matter in terms 
of an effective field theory with meson-nucleon coupling 
constants depending on the nucleon density and perform 
relativistic mean field or relativistic Hartree-Fock calcu¬ 
lations for finite nuclei using the semi-phenomenological 
density functionals, which were adjusted to reproduce 
DBHF results for infinite matter [sO, HH. For a recent 
discussion of such relativistic density functionals see [2^ 
and references listed there. Such studies consider the 
relativistic features of DBHF explicitly, but ignore more 
or less the effects of correlations beyond the mean field 
approach. It has been pointed out that the rearrange¬ 
ment terms originating from the density dependence of 
the meson-nucleon coupling constants are very important 
to improve the description of finite nuclei [3l|. 

In the present investigation we are going to discuss a 
different approximation scheme for the DBHF approach 
in finite nuclei. The focus is to consider the correlation 
effects in terms of a direct evaluation for finite nuclei 
and treat the Dirac effects in an approximate way. How¬ 
ever, in contrast to the method developed in [32| the 
Dirac effects are not described in terms of a local density 
approximation but are simulated using a three-nucleon 
interaction. 

Various attempts have been made to parameterize the 
3N forces discussed so far in terms of a simple local 3N 
force. As an example we mention the Urbana force [sM 
Is^l , which is composed of two terms 


V,,k = AV;^l + UV,, 


ijk 


''ijk ■ 


( 1 ) 


The first part is from 27r exchange with an intermedi¬ 
ate A excitation and may be considered to simulate the 
medium-dependent isobar effects discussed above. The 


second term is typically defined in terms of 2cr exchange 
and can be interpreted to simulate the effects of the Z- 
diagram discussed above. This means the second term is 
thought to represent the relativistic effects of the DBHF 
approach [ 3 ^. Typically this 3N force is reduced to a 
density-dependent NN interaction, which is then added 
to the bare NN interaction (see e.g. [s^ and references 
therein) and the parameters A and U in eq.(IT]) can be 
adjusted to reproduce the empirical saturation point for 
symmetric nuclear matter. 

This scheme has been criticized bey Hebeler and 
Schwenk[^ and later by Carbonne et al.jl^ They ar¬ 
gue that an expression for the total energy with kinetic 
energy ti, 2N interaction Vij and 3N potential Vijk 

E — ^ ) tiPi + — 'y ( VijPiPj -f — y ( VikjPiPkPj (2) 

i i,j 

leads to the single-particle energy 

-t- y ( ^ijPj T 2 'y ^ ^ikj PkPj : (H) 

3 3,k 

which is different from the result, which is obtained, when 
the 3N force is added to the 2N interaction by 

^Kf{p) = lv,, + lY.^.kjPk (4) 

k 

This is of course true and at first sight this would im¬ 
ply that the medium effects discussed above would lead 
to different results when they are treated in terms of a 
3-body force or considered as a density-dependent NN in¬ 
teraction. We note, however, that both approaches lead 
to the same result, if the single-particle energies are de¬ 
fined according to the Landau definition of the quasipar¬ 
ticle energy, i.e. 

= -^E{p ), (5) 

9 Pi 

which means that rearrangement terms due to the density 
dependence of are taken into account. With this 
inclusion the result is the independent on the treatment 
as a 3N term or a density-dependent 2N contribution. 

In this investigation we will discuss Brueckner-Hartree- 
Fock kind of calculations for nuclear matter and finite 
nuclei based on a realistic NN interaction with inclusion 
of a 3-nucleon force. Special attention will be paid to the 
rearrangement terms originating from the density depen¬ 
dence of the G-matrix and the treatment of the 3N force 
in terms of a density-dependent 2N interaction. We will 
show that an adjustment of the constant U defining the 
relativistic 3-body force in (HD is sufficient to obtain the 
empirical saturation point for symmetric nuclear matter. 
The same 3N force leads to a fair description also for the 
bulk properties of finite nuclei. 

After this introduction we will discuss the Brueckner 
Hartree Fock approach with 3N forces and inclusion of 
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rearrangement terms in section 2 of this paper. The re¬ 
sults for infinite matter and finite nuclei are discussed in 
section 3. Special attention will be paid to the treatment 
of the 3N force in finite nuclei and the description of the 
particle state spectrum in the Bethe-Goldstone equation. 
The main results and conclusions are summarized in sec¬ 
tion 4. 


II. BRUECKNER HARTREE EOCK AND 
REARRANGEMENT TERMS 

A. Nuclear Matter and Einite Nuclei 

The Brueckner Hartree Fock (BHF) approach can be 
defined in terms of three central equations. The first one 
of these equations is the Bethe-Goldstone equation 

G{uj) = V + V ^. G (6) 

Lu - Ho 

defining the so-called G matrix in terms of the free-space 
NN interaction V. Replacing the Pauli Operator Q, 
which forbids the scattering of the interacting nucleons 
into states, which are below the Fermi energy and there¬ 
fore occupied by other nucleons by the unit operator and 
the energy denominator w — Hq by the difference of ki¬ 
netic energies of free nucleons, the G matrix turns over 
into the Lippman Schwinger equation defining the scat¬ 
tering matrix T for two nucleons in the vacuum. There¬ 
fore the Bethe-Goldstone equation can be interpreted as 
the solution of the problem of two nucleons interacting 
in the nuclear medium and the G matrix can be under¬ 
stood as the effective interaction of two nucleons, which 
accounts for correlation between the interacting nucle¬ 
ons. The single-particle energies are then defined within 
the standard BHF using the Hartree-Fock expression in 
terms of the G-matrix interaction 

= {i\i\i) + ^(fj|G(w = ei+ ej)\ij)pj. (7) 

j 

The single-particle density pj is diagonal in the basis of 
Hartree-Fock states 


{k\p\j) = PjSjk , 


and the diagonal elements pj take the values 1 for occu¬ 
pied hole states with energies Sj below the Fermi energy 
and 0 for the particle states above the Fermi level. In 
the case of finite nuclei these Hartree Fock states have 
to be determined as the eigenstates of the BHF single¬ 
particle Hamiltonian, which corresponds to the single¬ 
particle energies defined in eq. 0. In the case of infi¬ 
nite nuclear matter these single-particle states are plane 
waves due to the symmetry of the system under trans¬ 
lational transformation. Note, however that also in the 
case of infinite matter a self-consistent solution of the 
Bethe-Goldstone equation (jS]) and the evaluation of the 


single-particle energies according eq. © is required to 
determine the starting energy w according to the Bethe- 
Brandow-Petchek theorem. 

After solving eqs.® and ([7]) in a self-consistent way 
one can evaluate the total energy as 

E = Y^{i\i\^)p^ + ^ Y^{ij\G{uj =e^+ ej)\ij)pjpi , ( 8 ) 

i ij 

which is the third of the 3 equations to define the BHF 
approximation. Note that using the BHF definition of the 
single-particle energies © the corresponding expression 
for the total energy leads to Koltuns sum rulef^ 


The definition of the single-particle energies for the 
intermediate particle states, i.e. the definition of the 
operator Hq in the propagator of the Bethe-Goldstone 
eq.® has been discussed for many years. It has been 
shown by Song et al. [d^ that the contribution from three- 
body correlations is minimized in nuclear matter with the 
so-called continuous prescription (4ll |. which means that 
the single-particle energies for the states above the Fermi 
level are calculated in the same way as those for the hole 
states below. As will be discussed below, we will try to 
adopt this prescription in our calculations. Since it is 
very elaborate to evaluate the single-particle energies for 
all states in finite nuclei, we will approximate the Hamil¬ 
tonian for the two particle states in this case by 

Ho = Q{ii+i2)Q-2G, ( 10 ) 

which is the operator of the kinetic energy of the inter¬ 
acting particles restricted to the states above the Fermi 
energy. A constant G is introduced to make the spec¬ 
trum “continuous” across the Fermi energy. Appropriate 
values will be discussed below. 

Another rather technical obstacle for the solution of 
the Bethe-Goldstone eq.® is the definition of the Pauli 
operator Q. This Pauli operator is easily defined for nu¬ 
clear matter using the rest-frame of the nuclear matter 
system by 


Q|fci, fe) 


|fci, fc 2 ) for |A:i| > kp and \k 2 \ > kp 
0 else. 


with ki denoting the momenta of the interacting nucle¬ 
ons and the Fermi momentum kp. The Bethe-Goldstone 
equation, however, is more easily solved in the center of 
mass frame of the interacting nucleon, as the momentum 
of the center of mass is conserved and the relative mo¬ 
mentum can be expanded in a partial wave basis. There¬ 
fore one typically employs the so-called angle-average ap¬ 
proximation for the Pauli operator and approximates the 
single-particle spectrum by a quadratic form 


e(fc) 


fc" 

2m* 


+ G* 


( 11 ) 
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with an effective mass m* and a constant C* fitted to 
describe the single-particle spectrum for the states be¬ 
low the Fermi momentum. Using these approximations, 
the Bethe-Goldstone equation can be solved separately 
in each partial wave, which reduces the numerical effort 
drastically. Methods have been developed to treat the 
Pauli operator and the single-particle spectrum without 
these approximations [421 l43l | and it has been shown that 
the angle-average in the Pauli operator is a reasonable 
approximation, while the parameterization of the single¬ 
particle spectrum according to eq. (HU can lead to con¬ 
siderable differences, as we will also discuss below. 

The problem of a precise treatment of the two-particle 
propagator is even more pronounced in calculation of 
finite nuclei, as the single-particle states are only de¬ 
fined after the corresponding Hartree-Fock equations 
have been solved. Since we are using the simple pa¬ 
rameterization of eq. (HU for the single-particle spectrum, 
we avoid a precise treatment for the Pauli operator ( 
see e.g.fa) and use the so-called angle-average for finite 
nuclei[3; Eg for a basis of oscillator states, which is ap¬ 
propriate for the nucleus under consideration. 

The single-particle states \i,lj) states are expanded in 
the very same oscillator basis \n,lj)HO 

\hlj) =^cl^\n,lj)HO, ( 12 ) 


assuming spherical symmetry of the states for the closed- 
shell nuclei considered. The quantum numbers I and j 
refer to the orbital and total angular momenta of the 
single-particle states and the projection quantum num¬ 
bers for the angular momentum are dropped. An attempt 
is made to optimize the oscillator basis in the sense that 
the oscillator parameter is chosen such that the corre¬ 
sponding expansion coefficients are close to 1 for all 
occupied single-particle states. In this way it is typically 
sufficient to restrict the expansion in eq. m to radial 
quantum numbers 0 < n < 4. 

The expansion coefficients ^ can then be determined 
by solving the BHF equations, which are given in the 
oscillator representation by 


{{n',lj\i+U^^^\nJj)HO+} Cni = : ( 13 ) 

n 


with (n', lj\i\n^ lj)HO the matrix elements for the kinetic 
energy in the basis of oscillator states and the correspond¬ 
ing matrix elements for the BHF single-particle potential, 
which can be calculated as 




BHF\ 


. O')i?o = X! 


2J- 


kl' j' mm' J 


2j + l 


i(n', Ij, m,' l'j'\G\n, Ij, m, l'j')HoPk 


(14) 


with the anti-symmetrized matrix-elements of the G- 
matrix {n'l'j'\G\n,lj,m,l'j')jjQ in the basis of 
two-nucleon oscillator states coupled to total angular mo¬ 
mentum J. The non-linear equations (HU and (HI are 
solved in an iterative way to obtain self-consistent solu¬ 
tions for the expansion coefficients as well as single¬ 
particle energies Si of the BHF single-particle states. 

The effects of the 3N interaction are taken into ac¬ 
count using a density-dependent 2N interaction as it is 
indicated in eq.(|ll) Note, however, that the weighting co¬ 
efficients of this density-dependent 2N interaction have 
been adjusted to obtain the correct expressions for the 
total energy and single-particle energies as presented in 
eqs. du and respectively. 

B. Rearrangement Terms 

The BHF approximation, which has briefly been 
sketched in the preceding subsection, corresponds to a 
Hartree-Fock calculation, replacing the two-particle in¬ 
teraction by the corresponding G-matrix. The G-matrix, 
however, must be understood as an effective interaction. 


due to its dependence on the starting energy uj and the 
Pauli operator, depends on the density operator of the 
system considered. Therefore the BHF definition of the 
single-particle energy does not obey the Landau defini¬ 
tion of the quasiparticle energy in ([5]). In fact, applying 
the Landau prescription to the energy functional (|5]) one 
obtains the BHF terms of © plus two additional terms, 
the starting energy rearrangement term A[/“ and the 
Pauli rearrangement term which are due to the 

dependence of G on starting energy ui and Pauli opera¬ 
tor Q. 

The starting energy rearrangement term can be writ¬ 
ten 

At/r = 

= '^PjPk{j,k\ — \j,k){i,j\G\i,j). (15) 

The second line of this equation is obtained by substitut¬ 
ing Ej in the hrst line by the corresponding BHF defini¬ 
tion of the single-particle energy. Note that adding AU^ 
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to the BHF definition of the single-particle energy leads 
to 


efBHF ^ + (16) 

= (* 1 ^ 1 *) + = Ei + : 

j 

(17) 

which means that we have replaced the single-particle 
density pj in eq.© by 


Pj — Pj 


1 + k) 


(18) 


This expression for Pj typically yields values of the order 
of 0.8... 0.9 and is often interpreted as a partial occu¬ 
pation of states j below the Fermi energy. The approxi¬ 
mation (IT6t represents the leading terms of the so-called 
Renormalized BHF approach (RBHF) ii- Therefore 
we will use this name also in the following. Note that 
the Koltun sum-rule of (0 cannot be used any longer to 
evaluate the total energy ([8]) using the RBHF definition 
of the single-particle energy. 

The Pauli rearrangement term can be written as 

= - E 10. • 

(19) 

and corresponds to the term of second order in G in the 
hole-line expansion of the self-energy. Calculations in¬ 
cluding Pauli- and starting energy rearrangement terms 
will be denoted as density-dependent Hartree-Fock cal¬ 
culations (DHF) and employ single-particle energies of 
the form 

e^HF ^ ^ 20 ) 


III. RESULTS 
A. Nuclear Matter 

All the calculations presented here have been per¬ 
formed using the proton-neutron part of the charge- 
dependent Bonn interaction (CD Bonn), which has been 
defined and adjusted to the 2-nucleon data by Machleidt, 
Sammarruca and SongQ. 

Results of conventional BHF calculations for symmet¬ 
ric nuclear matter, using the angle-average approxima¬ 
tion for the Pauli operator and the parameterization of 
the single-particle potential using the quadratic form of 
ea. ifTTI) are presented by the dashed line, labeled BHF 
AA, in Fig. [T] The single-particle potential, which is dis¬ 
played in the left panel of this figure for a Fermi momen¬ 
tum kp of 1.36 fm“^, which corresponds to the empiri¬ 
cal saturation density, reflects the quadratic parameter¬ 
ization, which is adjusted to reproduce the BHF single¬ 
particle potential U{k) for momenta k below the Fermi 
momentum and extended to momenta above kp. 


The calculated binding energy per nucleon of such BHF 
A A calculations, shown in the right panel of Fig. [1] yield a 
minimum, representing the prediction for the saturation 
point, at about twice the empirical saturation density 
and an energy of around -20 MeV, which is much more 
attractive than the empirical value of -16 MeV. 

The consistent treatment of the two-particle propaga¬ 
tor in the Bethe-Goldstone eq.®, avoiding the angle- 
average of the Pauli operator and using a consistent 
single-particle spectrum, leads to quite different results 
as can be observed from a comparison of the BHF results, 
presented by the red solid curves in Fig. [T] and the BHF 
AA results. As it has been discussed by Schiller et al.p^. 
these differences can mainly be attributed to the defini¬ 
tion of the single-particle potential. As can be seen from 
the left panel of Fig. [U the single-particle energies used 
to define the propagator of the Bethe-Goldstone equation 
are quite similar for momenta below kp. For the particle 
states with momenta above kp, however the calculated 
BHF energies are more attractive than described by the 
quadratic parameterization of the BHF AA approach. 

The corresponding differences in the two-particle prop¬ 
agator lead to matrix elements of G, which are in general 
more attractive in the BHF than in the BHF AA ap¬ 
proach, which leads to more binding energy in the former 
as compared to the latter calculation. This can be seen 
from the energy as a function of density curves, presented 
in the right panel of Fig. [T] The BHF calculations yield 
a saturation point with even larger binding energy (-24 
MeV) than the BHF AA approach at a larger saturation 
density. 

Fig. m also provides the relevant information about the 
effects of the rearrangement terms in the definition of the 
single-particle potential in eqs. (HU) and (1201) . The dom¬ 
inant contribution arises from the starting energy rear¬ 
rangement term, A{7“, which is taken into account using 
the RBHF approximation. As to be expected from the 
representation of the RBHF energies in terms of the par¬ 
tial occupation probabilities defined in dT51) these are less 
bound than the corresponding BHF single-particle ener¬ 
gies. This effect is more pronounced for the states with 
momenta below kp than for the particle states above the 
Fermi momentum. This enhances the calculated bind¬ 
ing energy even more, leading to a saturation point with 
an energy per nucleon of less than -25 MeV at a Fermi 
momentum of kp of 1.8 fm“^, which corresponds to a 
saturation density 2.3 times the empirical value. 

The Pauli rearrangement term, which is also included 
in the DHF approximation, has a small effect. Its main 
effect in the single-particle potential is concentrated at 
momenta around the Fermi momentum, where it leads to 
a reduction of the momentum dependence of U(k). This 
corresponds to the enhancement of the effective mass m* 
to the bare mass m, as it has already been discussed e.g. 
by Mahaux and Sartor and in [^. The weak contri¬ 
bution of the Pauli rearrangement effect is also reflected 
in the small difference of the energies calculated in the 
DHF as compared to RBHF approximation. 
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FIG. 1: (Color online) Results for symmetric nuclear matter calculations using the pn interaction of the CD Bonn potential. 
The left panel presents results for the single-particle potential U{k) assuming a Fermi momentnm kp of 1.36 fm“^, which 
corresponds to the empirical satnration densities. Resnlts are displayed for the BHF approximation nsing the angle-average 
in the Bethe-Goldstone equation (BHF AA) and BHF, RBHF (see eg. 11611 ') and DHF (see eg. 1201) ') calcnlations solving the 
Bethe-Goldstone equation without this approximation. The right panel shows the corresponding resnlts for the energy per 
nncleon calculated at various Fermi momenta. 


All the calculations discussed so far have been per¬ 
formed assuming just a realistic two-nucleon interaction, 
only, here the CD Bonn interaction, and we find, that 
the results of the calculated saturation points are part 
of the so-called Coester band[^ [T^. This is true for the 
BHF approach and the rearrangement terms just pro¬ 
vide a shift along this Coester band. As it is one of 
the main goals of this investigation, to simulate the rela¬ 
tivistic effects in terms of a 3N potential, we considered 
the 3N part of the Urbana interaction mo del [1^ [sj fix 
the parameter A in eq. (HD to be equal zero and adjust 
the parameter U, the strength parameter for the term to 
simulate the change of the Dirac spinors in the medium, 
to reproduce the empirical saturation point. We did not 
aim at a high-precision fit, but just tried to obtain results 
close to the experimental data. 

Results of such calculations are displayed in Fig. [5J It 
is worth noting that we can obtain a good description 
of the empirical saturation point by adjusting only one 
parameter, whereas most of the other attempts employ 
consider a 3 nucleon interaction with 2 or more parame¬ 


ters for a corresponding fit. 

The 3N Term (see eq.ljT])), which we consider in 
our studies, is of shorter range than the corresponding 
27r exchange term This may be the reason that 

the 3N term essentially provides a repulsive shift in the 
single-particle potential U{k) with almost no momentum 
dependence (see the example in the left panel of Fig. (2). 
So we do not obtain the strong momentum dependence 
which has been observed by Zuo et al. . 


B. Finite Nuclei 

The focus of the present investigation is to see if a 
BHF calculation with a parameterization of the Dirac 
effects in terms of a 3N interaction can provide a good 
description for the saturation point of nuclear matter as 
well as the bulk properties of finite nuclei. For that pur¬ 
pose we performed BHF and RBHF calculations for the 
closed shell nuclei and ^°Ca using the same 2N and 
3N interactions as just described for nuclear matter. 
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FIG. 2: (Color online) Results for symmetric nuclear matter calculations using the CD Bonn potential with (red curves) and 
without inclusion (black curves) of the 3N potential. The left panel presents results for the single-particle potential U{k) 
assuming a Fermi momentum kp of 1.36 fm“^ RBHF approximation. The right panel shows results for the energy per nucleon 
calculated at various Fermi momenta for the of BHF, RBHF and DHF approximations. 


Results of BHF calculations of using just the CD- 
Bonn potential are presented in table ID One finds that 
the calculated energy per nucleon (-6.08 MeV) is less at¬ 
tractive as compared to the experimental value (-7.98 
MeV) and the calculated radius for the charge distribu¬ 
tion, Rc, is much lower (2.35 fm) than the empirical value 
of 2.74 fm. In order to visualize this result for the “satu¬ 
ration point” for 160 in a way, which corresponds to the 
plot for nuclear matter as given e.g. in the right panel 
of Fig 121 we indicated this result in the energy versus the 
inverse of the radius of the charge distribution by a red 
dot in Fig. [H In fact, it is the upper of the two red dots, 
connected by a solid line in the left panel of this figure. 

Compared to the empirical data, represented by a 
green diamond, we observe a situation, which is quite 
different for ^®0 than for nuclear matter (see Fig. [5]). 
In both cases the BHF calculations yield value for kp or 
l/i?c, which are too large as compared to experiment, 
which implies that the average density calculated for the 
nuclear systems is too large. With respect to the energy 
per nucleon, however, the BHF calculations provide too 
much energy in nuclear matter and too little for the finite 
nucleus. Therefore in order to improve the comparison 
with experiment, the inclusion of the same 3N force must 


provide attraction in finite nuclei and repulsion in infi¬ 
nite matter and reduce the calculated saturation density 
in both cases. 

The inclusion of rearrangement terms going from the 
BHF to the RBHF approach yields occupation proba¬ 
bilities Pi of the order of 0.8 to 0.9 as shown in table 
m From eg. (TO it is obvious that this leads to less at¬ 
tractive single-particle energies in RBHF as compared to 
the BHF approach. The smaller attraction of the single¬ 
particle potential is reflected in a larger radius of the 
charge distribution. On the other hand the less attractive 
single-particle energies yield less attractive starting ener¬ 
gies u! in the Bethe-Goldstone equation, which leads to 
more attractive matrix elements of the G-matrix and re¬ 
sults in a more attractive energy per nucleon. This means 
that the inclusion of rearrangement terms shifts both, the 
energy per nucleon and the radius of the charge distribu¬ 
tion closer to the experiment. As one can see from table 
Uand Fig. [3]this effect is too small to provide a satisfying 
agreement. 

The study of nuclear matter discussed above already 
showed that the results of BHF kind of calculations are 
rather sensitive to a consistent treatment of the two- 
particle propagator in the Bethe-Goldstone equation. 
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NN only 

BHF RBHF 

with 3N 

BHF RBHF 

with 3N GD 

BHF RBHF 

Exp. 


e [MeV] 

e [MeV] 

P 

e [MeV] 

£ [MeV] 

P 

£ [MeV] 

£ [MeV] 

P 

e 


Protons 










Sl/2 

-58.19 

-48.69 

0.892 

-44.76 

-36.88 

0.917 

-41.59 

-32.89 

0.903 

-44 ± 7 

P3/2 

-27.05 

-20.93 

0.897 

-20.22 

-14.82 

0.840 

-17.64 

-12.11 

0.797 

-18.45 

Pl/2 

-20.02 

-16.25 

0.871 

-16.50 

-12.27 

0.824 

-14.20 

-9.80 

0.792 

-12.12 


Neutrons 










Sl/2 

-62.22 

-52.07 

0.892 

-48.36 

-39.67 

0.918 

-45.12 

-35.58 

0.907 

-47 

P3/2 

-30.98 

-24.19 

0.901 

-23.71 

-17.60 

0.846 

-21.02 

-14.69 

0.802 

-21.84 

Pl/2 

-23.83 

-19.44 

0.875 

-20.00 

-15.05 

0.829 

-17.61 

-12.36 

0.795 

-15.66 

E/A [MeV] 

-6.08 

-6.57 


-4.61 

-5.22 


-3.93 

-4.81 


-7.98 

Rc [fm] 

2.35 

2.45 


2.59 

2.66 


2.64 

2.72 


2.74 


TABLE I: Results for '^®0 using BHF and RBHF approximation without (NN only), with inclusion of the 3N interaction (with 
3N) and with 3N interaction in the global density approximation (with 3N GD). Values of single-particle energies (e) occupation 
probabilities {P, see ea llBII are listed for the occupied states as well as the energy per nucleon (E/A) and the radius of the 
charge distribution (Re)- The Pauli operator in the Bethe Goldstone equation has been defined in terms of oscillator function 
using an oscillator parameter b — 1.767 fm and a C = 5 MeV has been used to define the single-particle energies in eq. (nnj. 


The treatment of the particle-state spectrum in particu¬ 
lar requires special attention. The same is of course also 
to be expected for finite nuclei. As discussed in the pre¬ 
vious section, we approximate the Pauli operator Q in 
the Bethe-Goldstone eq.dH]) by a corresponding operator 
assuming oscillator states with an oscillator parameter b 
of 1.76 fm and & = 2 fm in the case of and ^°Ca, re¬ 
spectively. The spectrum of the particle states, denoted 
by Hq in eq.® is represented by the simple parameter¬ 
ization of (nni), trying to adjust the constant C in such 
a way that the resulting spectrum roughly matches the 
calculated single-particle energies. All the results dis¬ 
played in table H] have been evaluated with (7 = 5 MeV. 
A larger value of C leads to more attraction and so, the 
lines displayed in Fig. [3] indicate the range of results 
changing C from 5 to 15 MeV, with the less attractive 
energy representing the (7 = 5 MeV result. In this fig¬ 
ure the results of BHF calculations are represented by 
symbols connected by a solid red line, while the symbols 
connected by a dashed blue line refer to the results of 
RBHF calculations. The remarkable sensitivity of the 
calculation on this changes in the spectrum of the par¬ 
ticle states calls for more sophisticated investigations on 
this issue, to obtain unambiguous results. 

The effects of including the 3N force can be seen by 
comparing the results visualized in in Fig. |3] in terms 
of the square boxes with the corresponding results dis¬ 
played by circles. As to be expected, the inclusion of 
the 3N yields a repulsive effect and leads to a reduc¬ 
tion of the total energy accompanied by an increase of 
the nuclear radius. Comparing this effect with the corre¬ 
sponding repulsive effect one can obtain with a lowering 
of the parameter C just discussed, one finds that the 3N 
force yields a larger increase in the radius, if the energy 
is changed by a similar amount. This may reflect the fact 
that a proper treatment of the 3N force as compared to 


a modification of the 2N interaction has a larger effect 
on the single-particle potential (see the discussion in the 
introduction connected to eqs. and Therefore a 
repulsive 3N interaction will be more efficient in chang¬ 
ing the single-particle potential and the resulting radius 
of the particle distribution than a similar change in the 
2N interaction. 

Fig. [3] also shows results displayed in terms of crosses 
and denoted as 3N GD (see also table IJ) . As discussed 
above, this approximation scheme describes an attempt 
to evaluate the effects of the 3N force in nuclear mat¬ 
ter at various densities and then transfer the result to 
the calculation of finite nuclei by choosing the density 
of nuclear matter to be identical to the average particle 
density calculated for the nucleus considered. We find 
that this global density (GD) approximation yields ef¬ 
fects which are similar with a tendency to overestimate 
the corresponding effects of a direct treatment in the fi¬ 
nite nucleus by 20 to 30 percent. 

Similar calculations have also been done for the nu¬ 
cleus ^°Ca. Results on the “saturation properties” of 
this nucleus are displayed in the right panel of Fig. [3] 
considering the values (7 = 10 and (7 = 15 MeV for the 
parameterization of the particle state spectrum. More 
explicit results on the single-particle energies are shown 
in table El assuming (7 = 10 MeV. The main features of 
these results for ^°Ga are very similar to those discussed 
for the example and therefore confirm these findings. 

It is worth noting that, assuming (7 = 15 MeV, and 
including the 3N force RBHF calculations of both nuclei 
yield results for energy and radius of the charge distribu¬ 
tion, which are in a good agreement with the experimen¬ 
tal data (see Fig. [3] and table IHII) . We do not intent to 
celebrate this as a success of the Dirac BHF approach, 
or BHF approach simulating the Dirac effects of DBHF 
in terms of a 3N force. We will keep in mind that the 
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16 


o 


40 


Ca 




l/RJfm"'] 1/R^[fm‘‘] 


FIG. 3: (Color online) Results for the energy per nucleon and the radius of the charge distribution (Re) for (left panel) and 
“^^Ca are presented in a energy versus 1/Rc plot to enable the comparison with the corresponding figures for nuclear matter 
(Fig. [T]and Fig. [2]). Results referring to BHF calculations with different values of C in eq. m defining the particle spectrum 
in the Bethe-Goldstone equation are connected by a red solid line, those of RBHF calculations are connected by a blue dashed 
line. Results of calculations with just the 2N interaction, with inclusion of the 3N term and including the 3N term via a global 
density approximation are visualized by a dot, a box and a cross, respectively. The experimental result is given in terms of a 
green diamond symbol. 


strength of the 3N force, which was motivated to simu¬ 
late the Dirac effect, has been adjusted to reproduce the 
saturation point of nuclear matter. With the 3N force ad¬ 
justed we made a reasonable but not uniquely justified 
choice for the description of the particle-state spectrum 
to end up with a good description of the bulk properties 
of finite nuclei as well. 

Nevertheless, this result shows that using an appro¬ 
priate 3N force to simulate the effects of Dirac spinors 
modified in the nuclear medium within the framework 
of non-relativistic BHF calculations, one may be able to 
describe the bulk properties of nuclear matter and finite 
nuclei based on a realistic NN interaction. 

The description of bulk properties (energy, radius of 
particle distribution, density) is an important but only 
one of features of nuclear structure, which one hopes to 
describe within the relativistic DBHF. Other important 
aspects one hopes to describe within a relativistic de¬ 
scription of nuclear systems are the energy dependence 
of the optical potential[5l| and the spin-orbit splitting of 
the single-particle energies [s^, which is enhanced due to 


the enhancement of the small component of Dirac spinors 
in the nuclear medium. 


Is this enhancement of the spin-orbit splitting, which is 
important to describe the strength of the spin-orbit term 
observed in the experiment, also simulated by the simple 
3N force of eq. (HD? Inspecting e.g. the single-particle en¬ 
ergies of the P 3/2 and pi /2 states listed in table H] we do 
not find an enhancement of the spin-orbit splitting with 
the 3N force included. In fact, the differences between 
these single-particle energies are always smaller with in¬ 
clusion of the 3N force. This is related to the fact that 
the 3N force yields larger values for the radii which re¬ 
duces the spacing between the single-particle states. But 
even, if one accounts for this size effect, the 3N force 
does not provide an enhancement of the spin-orbit split¬ 
ting. This could be achieved by introducing an appro¬ 
priate spin structure in the 3N force, a feature which 
presumably would lead to more parameters and spoil the 
simplicity of the present approach. 
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NN only 


with 3N 


with 3N GD 


BHF 

RBHF 

BHF 

RBHF 

BHF 

RBHF 


e [MeV] 

£ [MeV] 

P 

£ [MeV] 

£ [MeV] 

P 

£ [MeV] 

£ [MeV] 

P 


Protons 









OSl/2 

-94.27 

-76.04 

0.850 

-53.18 

-36.65 

0.855 

-45.35 

-31.34 

0.790 

P3/2 

-63.36 

-48.11 

0.869 

-34.63 

-21.93 

9.745 

-27.58 

-17.92 

0.692 

Pl/2 

-53.82 

-41.61 

0.863 

-30.53 

-19.82 

0.723 

-24.10 

-16.16 

0.676 

<^5/2 

-33.20 

-21.81 

0.816 

-16.61 

-8.47 

0.718 

-11.36 

-5.51 

0.724 

lSl/2 

-26.61 

-18.06 

0.739 

-12.88 

-7.15 

0.726 

.9.33 

-5.04 

0.743 

d-3/2 

-18.22 

-12.20 

0.736 

-10.66 

-5.29 

0.713 

-6.33 

-2.85 

0.731 


Neutrons 









Sl/2 

-104.3 

-83.71 

0.854 

-61.35 

-42.02 

0.868 

-53.30 

-36.39 

0.810 

P3/2 

-72.74 

-55.31 

0.869 

-42.44 

-27.02 

0.760 

-35.02 

-22.71 

0.704 

Pl/2 

-63.12 

-48.74 

0.865 

-38.23 

-24.84 

0.737 

-31.41 

-20.89 

0.687 

^5/2 

-42.16 

-28.58 

0.830 

-23.99 

-13.36 

0.720 

-18.38 

-10.13 

0.722 

lSl/2 

-34.93 

-24.43 

0.752 

-20.13 

-12.00 

0.723 

-16.26 

-9.63 

0.736 

^3/2 

-26.88 

-18.72 

0.749 

-17.84 

-10.09 

0.712 

-13.19 

-7.39 

0.725 

E/A [MeV] 

-9.02 

-10.65 


-5.31 

-6.78 


-3.95 

-5.91 


Rc [fm] 

2.62 

2.72 


3.11 

3.38 


3.21 

3.52 



TABLE II: Results for ^°Ca using BHF and RBHF approximation. The Pauli operator in the Bethe Goldstone equation has 
been defined in terms of oscillator function using an oscillator parameter b = 2.0 fm and a C = 10 MeV has been used to define 
the single-particle energies in eq. dlOj. Further details see caption of table U 



NN only 

with 3N 



BHF RBHF 

BHF RBHF 

Exp 

i®0 E/A [MeV] 

-6.96 -9.47 

-5.67 -8.15 

-7.98 

Rc [fm] 

2.31 2.54 

2.52 2.77 

2.70 

E/A [MeV] 

-9.54 -12.40 

-5.98 -8.26 

-8.55 

Rc [fm] 

2.62 2.76 

3.08 3.43 

3.48 


TABLE III: Results for the energy per nucleon E/A and the 
radius of the charge distribution Rc for and ^°Ca calcu¬ 
lated in BHF and RBHF approximation are compared to the 
experimental data [i^, [S^. In contrast to the calculations 
leading to the results in table U and table El a shift G = 15 
MeV has been used to define the single-particle energies in 
eg.dTOl). 


IV. CONCLUSION 

An attempt has been made to simulate the relativistic 
features of the Dirac-Brueckner-Hartree-Fock (DBHF) 
approach by adding effects of a simple 3 nucleon (3N) 
force to non-relativistic many-body calculations based 
on the BHF approach. One parameter, the strength of 
the 3N force, is adjusted to reproduce the empirical sat¬ 
uration point of infinite nuclear matter and than used 


without further modifications for the description of finite 
nuclei. Special attention is paid to the energy spectrum 
of particle states, which is used in the propagator of the 
Bethe-Goldstone equation, and the importance of rear¬ 
rangement terms due to the energy dependence of the ef¬ 
fective interaction. For the studies of finite nuclei in par¬ 
ticular more effort is needed to optimize the treatment of 
the propagator in the Bethe-Goldstone equation. Taking 
a reasonable choice for the particle state spectrum and in¬ 
cluding the effects of rearrangement terms the Renormal¬ 
ized Brueckner Hartree Fock calculations with 3N force 
can reproduce the empirical values for the energy and 
radius of charge distribution of nuclei like and ^°Ca. 
This attempt to simulate the effects of the relativistic 
features of DBHF, however, fails to reproduce other pre¬ 
dictions of the Dirac phenomenology like the strength of 
the spin-orbit term in the single-particle field. 
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